Washington,  DC  20375-5320 


Naval  Research  Laboratory 


NRL/MR/5770— 12-9382 


Heave  and  Pitch  Dynamics  of 
Shallow  Draft  Surface  Vehicles 


Roger  S.  Cortesi 

Electronic  Warfare  Modeling  and  Simulation  Branch 
Tactical  Electronic  Warfare  Division 

Owen  Thorp 
George  Piper 

Systems  Engineering  Department 
United  States  Naval  Academy 
Annapolis,  Maryland 


April  10,2012 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
10-04-2012 

4.  TITLE  AND  SUBTITLE 


2. 


REPORT  TYPE 

Memorandum 


Heave  and  Pitch  Dynamics  of  Shallow  Draft  Surface  Vehicles 


6.  AUTHOR(S) 


Roger  S.  Cortesi,  Owen  Thorp,*  and  George  Piper* 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 

57-60 14-A95 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


NRL/MR/5770— 12-9382 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 


Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


NRL 


11 .  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

*Systems  Engineering  Department,  United  States  Naval  Academy,  Annapolis,  MD  21402-5002 


14.  ABSTRACT 

This  project’s  objective  is  to  develop  a  model  of  the  heave  and  pitch  dynamics  for  a  shallow  draft  surface  vehicle.  The  model  developed  uses 
a  novel  approach  of  representing  the  ship  dynamics  in  the  frequency  domain.  The  model  is  based  on  theoretical  first-principles  involving  highly 
nonlinear  equations  that  depend  on  the  water  wave  frequency.  Compared  to  time  domain  modeling  and  simulation  approaches,  the  frequency 
domain  model  provides  greater  insight  to  the  surface  vehicle’s  behavior  over  a  wider  operating  range.  To  validate  the  model,  experimental  data 
from  towing  tank  experiments  with  a  3  meter,  shallow  draft  boat  is  compared  to  theoretical  model  results.  The  experimental  data  compares  well 
with  the  model. 


15.  SUBJECT  TERMS 


Heave  Dynamics  Surface  vessel 

Pitch  Modeling  Wave-induced  motion 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Roger  Cortesi 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

UU 

24 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

(202)  767-6814 

1 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39. 18 


11 


Table  of  Contents 


Nomenclature . v 

1.0  Introduction . 1 

2.0  Review  of  Wave-Induced  Planar  Motion  Analysis . 1 

3.0  Systems  Analysis:  Frequency  Domain  Model . 4 

4.0  Experiment . 10 

5.0  Results . 13 

6.0  Conclusions . 17 

Acknowledgements . 17 

References . 18 


iv 


Nomenclature 

^  Free  surface  displacement  of  wave 

C,a  Wave  amplitude 


*4 

K 

V 

P 

z 

Za 

<k 


0 

0a 

0SS 

<t>e 

m 

az 

b 

c 

d 

e 

h 

lyy 


B 

C 

D 

E 

H 

8 

P 

K 

K 

B„ 


Wave  frequency 

Frequency  of  encountered  wave 

Wave  number 
Vehicle  velocity 

Relative  heading  of  vehicle  to  wave  train 
Heave  displacement 
Heave  amplitude 

Heave  phase  angle  relative  to  wave 
Heave  steady  state  displacement 
Pitch  displacement 
Pitch  amplitude 

Pitch  steady  state  displacement 

Pitch  phase  angle  relative  to  wave 
Vehicle  mass 

Added  mass  (the  mass  of  water  excited  by  the  vehicle's  motion) 

Heave  damping  coefficient 
Heave  stiffness  coefficient 
Pitch  inertia  coupling  term 
Pitch  damping  coupling  term 
Pitch  stiffness  coupling  term 
Pitch  axis  moment  of  inertia 

Added  mass  of  inertia 

Pitch  damping  coefficient 
Pitch  stiffness  coefficient 
Heave  inertia  coupling  term 
Heave  damping  coupling  term 
Heave  stiffness  coupling  term 
Acceleration  of  gravity 
Density  of  water 

Added  mass  coefficient 

Ratio  of  radiated  wave  amplitude  to  heave  motion  amplitude 
Vehicle's  sectional  beam  width 


v 


** ref 

an 

b 

n 

Cn 

F 

Fo 

Fs 

Fc 

/,./« 

M 

M0 

Ms 

Mc 

ms,  mc 

M 

C 

K 

x 

y 

u 

rt 


A  c 
I 

s 

G(s) 

Gzf(s) 

Gzm  (5) 

g8f  0) 
Gm  (^) 

GF(-(s) 


Sectional  reference  draft  of  the  wetted  surface 

Sectional  added  mass 

Sectional  damping  coefficient 

Sectional  stiffness  coefficient 

Wave  induced  force 

Amplitude  of  force 

Phase  angle  of  force  relative  to  wave 
Amplitude  of  in-phase  force  component 
Amplitude  of  out-of-phase  force  component 
Normalized  amplitudes  of  force  components 
Wave  induced  moment 
Amplitude  of  moment 

Phase  angle  of  moment  relative  to  wave 
Amplitude  of  in-phase  moment  component 
Amplitude  of  out-of-phase  moment  component 
Normalized  amplitudes  of  moment  components 

Inertia  matrix 
Damping  matrix 
Stiffness  matrix 
State  vector 
Output  vector 
Periodic  excitation  vector 
System  state  matrix 

System  input  matrix 

System  output  matrix 

Identity  matrix 
Laplace  transform  variable 

Transfer  function  matrix  of  heave  and  pitch  dynamics 
Transfer  function  from  force  to  heave 
Transfer  function  from  moment  to  heave 
Transfer  function  from  force  to  pitch 
Transfer  function  from  moment  to  pitch 
Transfer  function  from  wave  to  force 


VI 


Gzc(s) 

G6C(s) 

\Gz,  (7^e)  | 

^GZ(r(ja>e) 

\g0C(M)\ 

^G0C(jcoe) 


Transfer  function  from  wave  to  moment 

Transfer  function  from  wave  to  heave 

Transfer  function  from  wave  to  pitch 

Magnitude  of  transfer  function  Gz^(s)  at  frequency  COe 

Phase  of  transfer  function  Gz^(s)  at  frequency  COe 

Magnitude  of  transfer  function  G^(s)  at  frequency  COe 

Phase  of  transfer  function  G0^(s)  at  frequency  COe 


VII 


1.0  Introduction 

The  motivation  for  this  research  is  to  enhance  a  first  principles,  physics  based,  theoretical  model  to  replicate  small 
planing  boat  motions  as  described  in  Cortesi,  et.  al.  (2008)  and  Cortesi  and  Justh  (2007).  This  model  is  based  on 
semi-empirical  studies  of  planning  hulls,  such  as  those  developed  by  Savitsky  and  Gore  (1979)  as  well  as  the 
principles  discussed  in  Fossen  (1994)  regarding  the  modeling  of  displacement  hulls.  There  has  been  some  success 
in  using  this  relatively  simple  model;  but  it  was  deemed  desirable  to  enhance  the  features  of  the  model  by  adding 
the  effects  of  periodic  wave  action  on  the  small  boat.  Some  initial  success  was  achieved  in  this  endeavor  by  using 
strip  theory,  as  described  in  McCormick  (1973)  and  Bhattacharyya  (1978).  However,  the  first-principles  model  only 
worked  for  a  limited  range  of  boat  speeds  and  wave  frequencies.  It  was  soon  realized  that  for  all  of  its  simplicity, 
the  first  model  was  unable  to  incorporate  the  complex  effects  of  wave  action.  As  a  result,  the  NRL  model  was 
supplemented  for  one  which  allows  incorporation  of  the  non-linear  effects  of  a  small  boat  undergoing  motions  due 
to  wave  action  as  discussed  in  McCormick  (1973),  Bhattacharyya  (1978),  and  McCormick  (2010).  Additionally,  an 
approach  was  needed  which  allows  analytical  data  to  be  compared  to  experimental  data  over  the  full  range  of 

encountering  frequencies,  COe  .  It  was  therefore  decided  to  use  a  frequency  based  approach.  Using  this  approach, 
the  validity  of  the  analytical  model  was  able  to  be  determined  quickly  and  effectively. 

This  paper  is  organized  into  six  sections.  In  Section  1.0,  a  brief  introduction  is  given.  Section  2.0  presents  the 
fundamentals  of  wave  induced  planar  motions.  Section  3.0  discusses  the  coupled  equations  of  motion  for  heave 
and  pitch,  and  discusses  the  methods  used  to  determine  the  coefficients  for  the  equations  of  motion  and  the 
derivation  of  the  necessary  transfer  functions.  It  goes  on  to  discuss  the  derivation  of  the  vehicle  parameters  and 
the  derivation  of  the  wave  induced  force  and  moment  parameters.  Section  4.0  discusses  the  methodology  used  to 
conduct  the  validation  experiment.  Section  5.0  contains  the  results  of  the  model  and  the  experiment.  In  Section 
6.0,  a  summary  is  given  as  well  as  recommendations  for  future  work. 


2.0  Review  of  Wave-Induced  Planar  Motion  Analysis 

For  this  work,  it  is  assumed  that  the  vehicle  is  heading  into  the  waves  in  a  direction  that  is  transverse  to  the  wave 
crest  line.  The  seaway  consists  of  regular  harmonic  waves,  and  the  lateral  vehicle  motions  are  neglected.  Heave 
and  pitch  motions  are  referenced  to  the  vehicle's  body  axis  that  is  rigidly  located  at  the  vehicle's  center  of  gravity, 
CG.  The  planar  geometry  of  the  vehicle  is  illustrated  in  Figure  (1).  The  heaving  displacement,  z,  is  the  vertical 
displacement  of  the  CG  from  the  still  waterline;  while  pitch,  0,  is  the  angular  displacement  of  the  vehicle's  x-y 
plane  from  the  still  waterline.  The  free  surface  amplitude,  is  the  vertical  displacement  of  the  water  surface  from 
the  still  waterline  measured  at  the  vehicle's  CG. 


Manuscript  approved  January  5,  2012. 
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Vehicle  velocity 


The  effective  free-surface  displacement  encountered  by  the  vehicle  as  it  travels  is  described  by 

C  =  C  sin  (coet) 


(1) 


where  ^ a  is  the  wave  amplitude  and  COe  is  the  frequency  of  encounter.  This  frequency  is  related  to  the  vehicle's 
velocity,  V  ,  the  relative  heading  to  the  wave,  jd  ,  and  the  wave  frequency,  COw, according  to  the  relationship 


f 


co„  =  G)w 


cov 


1 - —  COS  jd 

V  8  J 


(2) 


For  the  case  where  the  vehicle  is  heading  into  the  wave  train,  the  encountering  angle  jd  is  180°. 


The  vehicle's  wave-induced  heave  and  pitch  motions  are  governed  by  the  vehicle's  inertial  interaction  with  the 
hydrostatic  and  hydrodynamic  forces  acting  on  it.  These  interactions  can  be  represented  by  two,  coupled,  second- 
order  differential  equations  with  periodic  forcing  functions,  as  discussed  in  McCormick  (1973),  Bhattacharyya 
(1978)  and  McCormick  (2010), 

( jyi  +  az  )z  +  bz  +  cz  +  dO  +  cO  +  hO  —  Fit )  (3.  a) 

(Iyy+Ayy)d  +  Bd  +  Cd  +  Dz  +  Ez  +  Hz  =  M(t)  (3.b) 


The  coefficients  in  Eqs.  (3. a)  and  (3.b)  correspond  to  the  inertia,  damping,  and  equivalent  stiffness  properties 
associated  with  the  vehicle-wave  interaction.  These  coefficients  are  a  function  of  the  vehicle's  particular  hull 
shape,  and  are  defined  in  Section  5.0.  The  terms  F(t )  and  M(t)  are  the  respective  wave-induced  excitation 
force  and  moment.  These  terms  will  be  discussed  Section  6.0.  Both  the  system  coefficients  and  the  excitation 
terms  are  dependent  on  the  encountered  wave  frequency,  COe .  In  addition,  the  excitation  terms  also  depend  on 

the  encountered  wave  amplitude,  ^ a .  Thus,  the  heave  and  pitch  dynamics  are  nonlinear  in  terms  of  frequency. 
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However,  if  we  constrain  our  discussion  to  a  particular  encountering  frequency,  CO e ,  of  interest  the  system 

coefficients  are  constant  and  Eqs.  (3. a)  and  (3.b)  can  be  considered  linear  with  regards  to  that  frequency. 
Therefore,  Eq.  (3)  can  be  conveniently  cast  in  compact  matrix  form,  as  in  Goldstein  (2001),  Meirovitch  (2010)  and 
Benaroya  (2004). 


M(oj) 


+  C(m) 


+  K(«) 


=  u (&>  ,0 


(4) 


where  the  matrices  are  defined  as  follows: 
The  inertia  matrix  is 

(m  + a)  d 


M(®> 

while  the  damping  matrix  is 

C(oj)  = 


D  (/w+^) 


b  e 
E  B 


and  the  stiffness  matrix  is 


K(tyJ  = 


c  h 
H  C 


Also,  in  Eq.  (4)  is  the  periodic  excitation  vector, 
U  (coe,t)  = 


F(coe,t ) 
M(a>  ,t) 


To  analyze  the  relationship  between  the  wave  induced  forces  and  the  heaving  and  pitching  motion,  the  two 
second-order  differential  equations  in  Eq.  (4)  are  formularized  into  a  fourth-order,  state  variable  model: 


x=rAx+rBu 


(5.a) 


and 


y  =  rcx 


(5.b) 


where  x  =  |^z  0  z  is  the  state  vector  representing  the  entire  state  of  the  system  at  any  given  time,  and 

y  =[z  #]  is  output  vector  denoting  the  variables  of  interest,  i.e.,  heave  and  pitch.  The  frequency  dependent 

systems  parameters  are  contained  in  the  system  matrices  defined  below. 
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0 

I 

0 

rA  = 

-M  K 

rB  = 

M  1 

3.0  Systems  Analysis:  Frequency  Domain  Model 

Knowing  that  the  heaving  and  pitching  dynamics  are  driven  by  the  periodic  wave  function  in  Eq.  (2),  it  is  highly 
advantageous  to  look  at  the  wave-induced  heave  and  pitch  motions  in  the  frequency  domain.  For  a  linear  system, 
such  as  the  one  represented  by  Eq.  (5),  a  sinusoidal  input  of  a  specific  frequency  results  in  an  output  that  is  also  a 
sinusoid  with  the  same  frequency,  but  with  a  different  amplitude  and  phase.  The  frequency-response  function 
describes  the  amplitude  change  and  phase  shift  as  a  function  of  frequency.  Thus,  for  a  given  seaway  characterized 

by  Eq.  (2)  with  a  wave  amplitude,  C,a  ,  and  a  frequency  of  encounter,  CO e ,  the  resulting  heaving  and  pitching 

motions  will  also  be  periodic  with  COe  and  characterized  by  amplitudes  Za ,  da  and  phase  angles  (j)z ,  (f)e  as  shown  in 
Figure  (2). 


Encountered  wave  motion  Resulting  vehicle  motion 

z  =  zasin(tf>/  +  ^;) 


*  0  -  0as\n{coet  +  (f>0) 


C  =  Ca  Sin(fi>eO 


Longitudinal 

Dynamics 


Fig.  2  Input/Output  Relationship  between  Wave  and  Vehicle  Motion 

To  obtain  a  frequency-response  function,  employ  the  transfer  function  concept  which  relates  the  system  output  to 
the  system  input  in  the  Laplace  domain.  The  transfer  function  concept  allows  one  to  convolve  input/output 
relationships  algebraically,  as  in  Meirovitch  (1985).  In  the  Laplace  domain,  Eq.  (5)  can  be  written  in  a  compact 
input-output  relation  as 


yO)  =  G(s)u(,s)  (6) 

where  y (s)  and  u(s)  are  the  Laplace  transforms  of  the  output  and  input  y(t)  and  u(^),  respectively.  The  term 
G^)  is  a  transfer  function  matrix,  which  is  an  algebraic  expression  defined  by  the  ratio 


G(s)  - 


y(s) 

00) 


(7) 


The  transfer  function  matrix  G(s)  can  be  obtained  from  Eq.  (5)  by 


G(s)  =  Tc[sl-TA]-lTB 

_  Gzf  C5)  GZm 

_G0f(s)  G6m(s) 


4 


The  elements  of  the  matrix  G(s)  are  scalar  transfer  functions  representing  the  interactions  between  heave  (z), 

z(s') 

pitch  [0),  the  wave  induced  force  (F)  and  moment  (M)  as  defined  by  GZF(s)  = - ,  GZM(s)  = - , 

F(s )  M(s ) 

^  ,  0(s)  „  .  .  0(s) 

GaF  (5)  = - ,  Gm  ( s )  = - ,  respectively.  Therefore,  the  heave  and  pitch  responses  represented  by  Eq. 

F(s)  M(s ) 

(6)  can  be  written  more  explicitly  as 

z(s)  =  Gzf(s)F(s)  +  Gzm(s)M(s )  (g) 

6(s)  =  G0F(s)F(s)  +  Gm(s)M(s )  (10) 

The  relationship  between  the  respective  wave-induced  force  and  moment  and  the  encountered  waves  can  also  be 
expressed  in  the  Laplace  domain  as 

F(s)  =  Gf((s)C(s)  (11) 

M(s)  =  Gm^s)C(s)  (12) 

^  F(s)  M(s ) 

where  transfer  functions  C Trr{s)  = -  and  GM/-{s)  = -  represent  the  dynamics  of  the  wave  force  and 

moment  induced  at  the  encountered  waves.  The  derivation  of  the  force  and  moment  transfer  functions,  GF ^  (s) 
and  Gmc  (s),  is  discussed  in  Section  6.0. 

By  substituting  Eqs.  (11)  and  (12)  into  Eqs.  (9)  and  (10),  the  heaving  and  pitching  motions  are  now  related  to  the 
encountered  waves  by 


zO)  =  \GZF  (s)Gfc  (s)  +  Gzm  (s)Gmc  0)]  £  0) 


m = [g6F  (. s)gf(  (s) + cm  {s)gm(  g)\ as) 


respectively. 

Expressing  these  relations  in  terms  of  transfer  functions,  one  obtains 


G7As)  = 


as) 

=  Gzf  (, s)Gf(  (5)  +  Gzm  (, s)Gm(  (5) 
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(16) 


Gg((s)  = 


0(s) 

as) 


=  GeF{s)GF((s)  +  GeM{s)GMas) 


OM 


For  a  sinusoidal  wave  input,  where  ^  =  £aSin ((Det) ,  it  can  be  shown  that  when  steady  state  conditions  are 
reached,  the  heaving  and  pitching  frequency  responses, 

zJt)  =  zasm(G)et  +  <f>z) 

ejt)  =  Oii  sin(®/  +  4  ), 


can  be  calculated  by  replacing  the  Laplace  variable  s  in  the  transfer  function  Eqs.  (15)  and  (16)  by  jo)e  such  that 


Za  =  \GZ(  ( jcoe )  </>z=  ZGz( (M ) 

(17) 

da  =|G^(M)|C'  4  =  zGe((jcoe) 

(18) 

Here,  j  =  V— 1  . 

On  this  basis,  and  are  amplitude  ratios  of  heave  and  pitch  resulting  from  the 

encountered  wave.  Similarly,  AGZ(,{jcoe)  and  AG9^{jcoe)  are  the  heave  and  pitch  phase-shift  with  respect  to 

the  encountered  wave.  The  amplitude  ratios  and  phase  shifts  are  functions  of  the  encountered  wave  frequency 
and  therefore,  can  be  conveniently  presented  in  graphical  form  as  illustrated  later  in  the  paper.  These  figures 
show  the  heave  amplitude  ratio  versus  frequency  of  encounter.  Given  the  wave  amplitude  and  frequency  for  a 
specific  seaway,  the  resulting  heave  amplitude  can  be  read  directly  from  such  a  graph. 

3.1  Vehicle  Parameters 

The  inertia,  damping,  and  stiffness  coefficients  in  Eq.  (3)  are  defined  in  terms  of  strip  theory  as  presented  by 
Krovin-Kroukovsky  and  Jacobs  (1957).  With  this  method,  the  vehicle  is  divided  into  a  finite  number  of  transverse 
sections  along  its  length.  The  inertia,  damping,  and  stiffness  coefficients  for  each  section  are  considered  constant 
and  are  denoted  by  the  subscript  n.  The  overall  vehicle  coefficients  are  determined  by  integrating  the  section 
coefficients  over  the  vehicle  length,  L.  The  definitions  of  Eq.  (3)  coefficients  are  summarized  in  Table  1. 

Table  1  Hydrodynamic  Coefficient  Definitions 


Term 

Definition 

m  =  j 

L 

m„d % 

Vehicle  mass 

az=. 

f  and% 

L 

Added  mass  (the  mass  of  water 
excited  by  the  vehicle's  motion), 

6 


L 

\d% 

Heave  damping  coefficient 

L 

'M 

Heave  stiffness  coefficient 

d  =  - 

{  an£d4 

L 

Pitch  inertia  coupling  term 

e  =  — 

bn%d%  +  uaz 

Pitch  damping  coupling  term 

h  =  - 

[  cn%d%  +  ub 

L 

Pitch  stiffness  coupling  term 

I„  = 

f  mj2d^ 

L 

Pitch  axis  moment  of  inertia 

K  = 

\aj2d^ 

L 

Added-mass  of  inertia 

L 

bj2d{ 

Pitch  damping  coefficient 

C=1 

L 

cJ2d?-uE 

Pitch  stiffness  coefficient 

D  =  d 

Heave  inertia  coupling  term 

E  =  - 

\bjd^-uaz 

L 

Heave  damping  coupling  term 

H=- 

-\CnZd% 

L 

Heave  Stiffness  coupling  term 

The  term  mn  in  Table  1  is  sectional  vehicle  mass.  The  term  an  is  the  sectional  added  mass  and  can  be 
determined  from 


an='-K,PXBl 


(19) 


where  p  is  the  mass  density,  Bn  is  the  vehicle's  sectional  beam  width,  and  ka  is  the  added  mass  coefficient.  ka 
is  a  function  draft/beam  ratio,  shape  of  the  hull's  cross  sectional  area,  and  the  wave  frequency,  .  Values  of  ka  for 

various  hall  cross  section  shapes  can  be  found  in  Lewis  (1929).  The  term  bn  is  the  sectional  damping  coefficient 
and  can  be  determined  from 


b  = 


pill 

col 


(20) 


where  g  is  the  gravitational  constant  and  Rz  is  the  ratio  of  radiated  wave  amplitude  to  heave  motion  amplitude. 
The  expression  for  Rv  used  in  this  work  was  derived  in  Yamamoto  et.  al.  (1980)  and  is  given  by 
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R=2e  8  sin 


r  B  co 2A 

n  e 

Vs) 


(21) 


The  term  zref  is  a  sectional  reference  draft  of  the  wetted  surface. 

The  sectional  stiffness  coefficient,  Cn  ,  is  associated  to  the  buoyancy  restoring  force  due  to  heave  disturbances  and 
is 


c„  =  PSBn 


(22) 


3.2  Wave  Induced  Force  and  Moment 

In  general,  the  force  induced  by  harmonic  waves  can  be  written  as 

F(t)  =  Fc  cos (coet)  +  Fs  sin(<»/) 

=  F0  sm(coj  +  0F) 


(23) 


where  F0  is  the  force  amplitude,  defined  by 


F0=s[fF^ 


(24) 


and  (j)F  is  the  phase  shift  of  the  force  with  respect  to  the  wave.  That  is, 


(f)F  =  tan 


-i 


rF^ 

vTy 


(25) 


The  force  amplitude  terms  Fc  and  Fs  are  derived  in  Goldstein  (2001),  Meirovitch  (2010)  and  Benaroya  (2004) 
and  can  be  expressed  as 


FM)  =  Cae~  zfM) 


fc  W  =  j ( cn  -  °fean ) sin (K^)d^  +  coej(bn-u  cos (K^)d^ 

L  L  dC, 


(26) 

(27) 


Fs(coe)  =  Cae-KzfM ) 


(28) 


fs (®e )  =  J ( cn  - ) cos (K^)d^  —  coej (bn  -u  ^f) sin (K^)d^ 

L  L 


(29) 


Substituting  Eqs.  (26)  and  (28)  into  Eq.  (24)  and  dividing  the  resulting  equation  by  the  wave  amplitude,  one  obtains 
an  expression  for  the  amplitude  ratios  of  force  to  wave.  That  is, 
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Y^(°e)  =  e~ 

*2  a 


Kz  r  2  ,  rl 


Similarly,  substituting  Eqs.  (26)  and  (28)  into  Eq.  (25)  one  obtains  an  expression  for  the  phase  shaft  between  the 
force  and  wave  which  is, 


=  tan”  -f 

V  j  s 


The  frequency  transfer  function  between  the  wave  induced  force  and  the  encountered  waves  can  now  be  written 


GFcUa>e)  =  -re 


-  - 0  jyj&F 


In  a  similar  fashion,  the  moment  induced  by  harmonic  waves  can  be  written  as 

M(t )  =  Mc  cos (coet)  +  Ms  sin(<»/) 

=  M0  sin (o)et  +  </)M) 

where  M0  is  the  moment  amplitude  defined  by 


Mo=4Mc+Ms 


id  (j)M  is  the  phase  shift  of  the  moment  with  respect  to  the  wave,  and  is  expressed  as 


The  moment  amplitude  terms  M  and  M  can  be  expressed  as 


Mc{oJe)  =  ^aeKzmc{o}e) 


mc  (»«)  =  /  (cn  -  )%  sin(X£)^  +  a>ef  (bn  -u-?f)%cos(K%)d% 

L  L 


MM)=£ae'KZmM) 


*s (®e )  =  J ( c„  - co2ean )4 cos (K^)d^  - coe  J (bn  - u sin (K^)d^ 

L  L 
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Substituting  Eqs.  (36)  and  (38)  into  Eq.  (34)  and  dividing  the  resulting  by  the  wave  amplitude,  an  expression  for  the 
amplitude  ratios  of  moment  to  wave  is  obtained.  That  expression  is 


^{coe)  =  e~Kz  <Jm[+m2s  (40) 

a 

Similarly,  substituting  Eqs.  (36)  and  (38)  into  Eq.  (35)  an  expression  for  the  phase  shaft  between  the  moment  and 
wave  is  found  to  be. 


<j)Jcoe)  =  Vcm 


\ms  / 


(41) 


The  frequency  transfer  function  between  the  wave  induced  moment  and  the  encountered  waves  can  now  be 
written  as  follows: 


Gu(Uo>,)  =  ^e> 

^  a 


(42) 


Eqs.  (32)  and  (42)  can  now  used  with  Eqs.  (16)  and  (17)  to  develop  the  frequency  response  model. 


4.0  Experiment 

To  validate  the  model  developed  above,  the  wave-structure  interaction  of  a  small,  shallow-draft,  flat-bottom  boat 
was  investigated  in  the  380  foot  towing  tank  of  the  U.  S.  Naval  Academy's  Hydromechanics  Laboratory,  described 
in  http://usna.edu/Hydromechanics/homepage.html.  The  test  boat  had  the  following  characteristics: 

Length  =  3.05m 

Beam  =  1.142  m  at  L/2 

Mass  =  49.9  kg 

To  simulate  the  effects  of  onboard  personnel  and  ancillary  equipment,  120  kg  of  weights  were  added  in  the  aft 
section  of  the  boat  as  sketched  in  Figure  3.  Additionally,  the  boat  was  attached  at  its  center  to  the  towing  carriage 
via  a  set  of  heave  posts  which  rode  freely  when  any  waves  were  encountered.  The  overall  mass  of  the  heave  post 
was  17.4  kg,  making  the  overall  mass  of  the  boat  187.3  kg.  Figure  4  shows  a  schematic  of  the  Naval  Academy 
Towing  Tank. 
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Max  Beam  at  Gunnels  =  44.8” 


Fig.  3  Schematic  of  Test  Boat 


Plan  View 


Towing  Carriages 

I - II - 1 


t 


Wave  maker 


Elevation  View 


Wave  Absorber  -  Beach 


Maximum  Speed:  32  ft/sec  (9.8  m/sec) 
Maximum  Wave  Heiaht:  1.5  ft  (0.46  m) 


Fig.  4  Schematic  Showing  Plan  and  Elevation  Views  of  Naval  Academy  Tow  Tank 

The  380-foot  towing  tank  is  large  enough  to  test  ship  models  up  to  25  feet  in  length  and  weighing  several  thousand 
pounds.  The  tank  length  allows  testing  at  high  speeds.  The  tank  is  outfitted  with  two  towing  carriages,  a  dual-flap 
wavemaker  and  specialized  equipment  for  measuring  resistance,  seakeeping  and  maneuvering  characteristics  of  all 
types  of  marine  vehicles  and  ocean  platforms.  The  test  boat  was  attached  to  the  towing  carriage  as  seen  below  in 
Figures  5  and  6. 
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Fig .  5  Photograph  of  the  Test  Boat  Connected  to  Towing  Carriage  in  Towing  Tank 


Fig .  6  View  of  the  Test  Boat  from  the  Carriage 

Twenty  experimental  runs  were  made,  while  varying  three  parameters:  wave  height,  wave  frequency,  and  ship 
speed.  A  summary  of  the  experimental  runs  is  given  in  Table  2. 
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Table  2  Experimental  Heave-to-Wave  Amplitude  Ratios 


Boat 

Velocity 

(Kts) 

Wave 

Frequency 

(Hz) 

Encountered  Wave 
Frequency 
(Hz) 

(za) 

Measured  Heave-to-Wave  Amplitude  Ratio,  — 

vC  7 

1"  Wave  Height 

2"  Wave  Height 

3"  Wave  Height 

3 

0.4 

0.56 

0.992 

1.088 

1.105 

3 

0.7 

1.18 

0.752 

0.712 

; 

3 

1.0 

1.99 

0.140 

0.109 

- 

4 

0.4 

0.61 

1.035 

1.085 

1.141 

4 

0.7 

1.35 

0.510 

0.723 

i 

4 

1.0 

2.32 

0.152 

0.180 

- 

6 

0.4 

0.72 

1.196 

1.361 

1.397 

6 

0.7 

1.67 

0.885 

0.860 

- 

6 

1.0 

2.98 

- 

0.275 

- 

5.0  Results 

In  this  section,  experimental  results  are  used  to  validate  the  dynamic  model  developed  herein.  The  boat  heave 
and  wave  height  for  each  experimental  run  was  measured  and  recorded  in  time.  To  illustrate  the  experimental 
results,  the  heave  time  responses  at  boat  speeds  of  3,  4,  and  6  Kts  to  1  inch,  0.4  Hz  harmonic  waves  are  shown  in 
Figures  7  through  9. 


0.05 


10  12 

Time  (sec) 


Fig .  7  Heave  Time  Response  to  1  Inch  -  0.4  Hz  Waves,  Vehicle  Velocity  =  3  Kts 
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Fig .  8  Heave  Time  Response  to  1  Inch  -  0.4  Hz  Waves,  Vehicle  Velocity  =  4  Kts 


Time  (sec) 


Fig.  9  Heave  Time  Response  to  1  Inch  -  0.4  Hz  Waves,  Vehicle  Velocity  =  6  Kts 
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The  experimental  results  show  that  the  heave  response  to  regular  harmonic  waves  is  also  harmonic  with  frequency 
equal  to  the  encountered  wave  as  expected.  Also  shown  in  Figures  7  through  9  are  results  of  heave  time  response 
simulations  obtained  from  the  heave  dynamic  model  in  Eq.  (13).  The  measured  time  history  of  the  wave  height 
was  used  as  input  for  the  simulations.  As  can  be  seen,  the  time  simulation  and  the  experimental  results  compare 
well. 

Knowing  that  the  boat's  response  to  regular  harmonic  waves  is  also  harmonic,  it  is  highly  advantageous  to 
represent  the  boat's  behavior  as  a  frequency  response.  The  frequency  response  conveniently  characterizes  the 
boat's  response  to  a  wider  range  of  seaway  conditions,  as  opposed  to  time  responses.  From  the  experimental 

data,  the  ratio  of  the  boat  heave  amplitude,  za ,  and  the  wave  amplitude,  <^a ,  were  calculated  and  shown  in  Table 
2. 

The  experimental  heave-to-wave  amplitude  ratios  are  plotted  as  a  function  of  the  frequency  of  encounter  ( COe ) 

for  the  three  boat  velocities  in  Figures  10  to  12.  The  figures  show  that  the  heave  amplitude  ratio  is  fairly 
independent  to  wave  height  changes  at  a  given  frequency.  This  result  supports  the  linearity  assumption  of  the 
dynamics  described  in  Eq.  (3)  at  a  specified  input  frequency. 


Fig.  10  Heave  Frequency  Response  to  Wave  Input,  Vehicle  Velocity  =  3  Kts 
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Fig.  11  Heave  Frequency  Response  to  Wave  Input,  Vehicle  Velocity  =  4  Kts 


Fig.  12  Heave  Frequency  Response  to  Wave  Input,  Vehicle  Velocity  =  6  Kts 
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The  frequency  response  amplitude  ratio  |Gz^(j^)|  in  Eq.  (17)  was  calculated  for  the  test  boat,  and  is  shown 

superimposed  on  the  experimental  data  in  Figures  10  through  12.  The  frequency  domain  model  shows  the 
characteristics  associated  with  a  fourth-order  linear  physical  system.  That  is,  constant  low  frequency  amplitude 
with  resonate  peaks  around  the  natural  frequency,  and  a  -40  dB  per  decade  slope  at  high  frequencies.  The  heave 
natural  frequency  can  be  calculated  by 


=  271  fn  =  . 


'm  +  a„ 


(43) 


For  the  test  boat,  the  natural  frequency  is  approximately  equal  to  fn  =  0.64  Hz.  The  peaks  and  valleys  observed  at 
frequencies  greater  than  0)n  correspond  to  nonlinearities  in  the  buoyancy  force  and  moment.  Buoyancy  depends 

heavily  on  the  relationship  between  wave  length  and  boat  length.  Buoyancy  will  decrease  when  the  boat  length  is 
an  integer  multiple  of  the  wave  length  and  will  increase  when  it  is  not.  Therefore,  as  the  wave  frequency  increases 
and  the  wave  length  decreases,  the  buoyancy  force  and  moment  will  rise  and  fall. 

In  general,  the  frequency  domain  model  corresponds  well  with  the  experimental  data,  especially  at  low 
frequencies.  Any  discrepancies  observed  between  the  model  and  data  can  be  attributed  to  inaccuracies  in  the 
model's  added  mass  and  damping  coefficient  terms.  Expressions  for  these  terms  were  derived  from  interpolations 
of  experimental  data  and  therefore,  inherently  have  a  degree  of  uncertainty.  In  addition,  experimental  error 
associated  with  the  towing  tank  experiments  may  have  also  contributed  to  any  discrepancies. 


6.0  Conclusions 

In  this  paper,  a  model  for  the  coupled  heave  and  pitch  dynamics  of  surface  vehicles  is  developed  using  a  novel 
approach  of  representing  the  ship  dynamics  in  the  frequency  domain.  The  model  is  based  on  theoretical  first 
principles  involving  highly  nonlinear  equations  that  depend  on  the  wave  frequency.  Compared  to  time  domain 
modeling  and  simulation  approaches,  the  frequency  domain  model  provides  greater  insight  to  the  surface  vehicle's 
behavior  over  a  wider  operating  range.  The  theoretical  model  was  validated  by  experimental  data  that  was 
obtained  from  towing  tank  experiment.  The  test  boat  has  a  3-meter  shallow  draft  and  a  flat  bottom.  The 
theoretical  model  results  compare  well  with  the  experimental  data. 
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